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ABSTRACT 

By explicitly taking into account effects of Alfven waves, I derive from a simple energetics argument 
a fundamental relation which predicts solar wind (SW) speeds in the vicinity of the earth from physical 
properties on the sun. Kojima et al. recently found from their observations that a ratio of surface 
magnetic field strength to an expansion factor of open magnetic flux tubes is a good indicator of 
the SW speed. I show by using the derived relation that this nice correlation is an evidence of the 
Alfven wave which accelerates SW in expanding flux tubes. The observations further require that 
fluctuation amplitudes of magnetic field lines at the surface should be almost universal in different 



o 

. coronal holes, which needs to be tested by future observations. 
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1. INTRODUCTION 

Speeds of the solar wind (SW) in the vicinity of the earth vary from -300 to -800km/s ()Phillips et al.lll995(l . The 
SW speed is one of the important parameters to predict geomagnetic storms triggered by interactions between the SW 
>>■ ' plasma and the earth magnetosphere (e.g. Wu & Lepping 2002). If we can tell SW conditions near the earth from 
, observed properties on the sun, we can forecast geomagnetic conditions beforehand since it takes a few days until the 
SW reaches us after emanating from the sun. 
^— ^ Thus, various attempts have been carried out to derive simple relations which connect physical properties on the 
^ I sun and SW speeds near the earth. Wang & Sheeley (1990; 1991, hereafter WS90 and WS91) showed that SW speeds 
are anti-correlated with expansion factors of magnetic flux tubes from their long-term observations as well as by a 
simple theoretical model, and this relation is widely used to predict SW speeds (e.g. Arge & Pizzo 2000). Fisk, 
Schwadron, & Zurbuchen (1999; here after FSZ) claimed that t he SW speed should have a positive dependence on the 
magnetic field strength on the sun. ISchwadron fc McComas I (|2003..) puts forward a SW scaling which explains the 
Oh| obser ved anti-correlati on between the SW speed and freezing-in temperatures, reflecting the coronal temperature, of 
Q . ions l)Geiss et al.lll995l) . iMcIntosh fc Leamon I l)2flfl,'Tf l further introduced a correlation between a scale length in the 
^ chromosphere and the SW speed. 

^ Turning to the acceleration mechanism of the SW, it is generally believed that the Alfven wave is a promising 
■ candidate which dominantly works both in heating and accelerating the SW plasma (Belcher 1971; Ofman 2004, 
Cranmer 2005; Suzuki & Inutsuka 2005; 2006; hereafter, SI05 and SI06). However, there is no fundamental relation 
derived so far, which is directly linked with the Alfven wave. The aim of the present paper is to derive a simple formula 
which connects the SW speed and the solar surface conditions through Alfven waves by referring to results of recent 
^ ' numerical simulations (SI05, SI06). 

^ From an observational viewpoint Kojima et al.(2005) have extensively surveyed relations between SW speeds around 
one astronomical unit (lAU), wiaUj and properties of magnetic flux tubes, radial magnetic field strength at the 
photosphere, -Br,©, and a (total) super-radial expansion factor of the tube, /tot during a solar minimum phase, 1995- 
1996 (Figure 1), where the values are averaged over each open coronal hole and the potential field-source surface method 
(e.g. Sakurai 1982) is used to derive /tot- They claimed that a ratio, i?r,o//tot, is the best indicator of uiau (top panel) 
(see also Suess et al.l984), whereas they also found a moderate correlation of wiau — l//tot (middle panel) (WS90) 
and a weak correlation of viav — -Br,© (bottom panel) (FSZ). Note that, only within the framework of the potential 
field-source surface method, -Br,©//tot is equivalent with magnetic field strength at the source surface (assu med at 
2.5-Rf 7i; Rq is the solar radius), the outside of which field lines are assumed be radially oriented (Hakamada & K oiimal 
[1999"), while they use the ratio of -Br,© and /tot as it stands because it is more physically motivated (see §3). /tot used 
in this letter is defined as the total expansion factor from the solar surface to lAU. 

The obtained nice correlation of wiau — -Br,©//tot seems quite reasonable; the positive correlation on -Br,© appears 
natural since -Br,© controls strength of Poynting energy which is injected from the surface and finally accelerates 
the SW (FSZ); the negative dependence on /tot seems reasonable as well because /tot determines adiabatic loss of 
the SW in the flux tubes (WS91) (Figure 2). One may further speculate that the wiau — Sr.©//tot relation reflects 
Alfven waves, a type of Poynting flux, in the diverging flux tubes. Here I develop this consideration to give a quantitate 
interpretation of the relation. 



Electronic address: |stakeru@scphy s . kyoto-u . ac. jp | 

^ Department of Physics, Kyoto University, Kitashirakawa-Oiwake-cho, Sakyo-ku, Kyoto, Japan 606-8502; JSPS Research Fellow 



2 



T. K. Suzuki 



CO 



900 


1 1 1 


800 




/UU 




600 




oUU 




400 






- /♦ ' 


300 




200 


1 i/ 1 




J I I I I I I I L. 







0.5 1 
B,0Aot(G) 



1.5 



900 
800 
700 
5 600 
< 500 
^ 400 
300 
200 



900 
800 
700 
^ 600 
< 500 



-| 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



4^* ♦ 



♦ V 



♦.• ♦ ♦♦ 







0.1 0.2 

Vftot 



0.3 



- m 



400 
300 
200 



-1 — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — I — r 



I* 



I— — ♦i=*^r* — ' 



I ♦ 1 



I — ♦ — I 







10 
B 



15 



20 



25 



r.o 



(G) 



Fig. 1. — Relations between SW speeds at r ~ lAU, viaV} and properties of magnetic flux tubes. Observed data are from Kojima 
et al . (2005). Coronal magnetic fields are extrapolated from -Br,0 by the potential field-source surface method (Hakamada & Kolirn^ 
119991) . /tot is derived from comparison between the areas of open coronal holes at the photosphere and at the source surface (r = 2.5-R0). 
■flAU is obtained by interplanetary scintillation measurements. diaUi ^r-,©, and /tot are averaged over the area of each coronal hole and 
the data points correspond to individual coronal holes. (Top) : diau on i?r,© //max- Lines are theoretical prediction from Equation Bl . 
Solid line indicates the fiducial case {{5B±5v±) = 8.3 X lO^G cm s~^ and Tq = 10®K). Dot-dashed line adopt higher coronal temperature 
(Tc = 1.5 X lO^K) with the fiducial {5B±Sv±). Dashed line adopt smaller (SB±Sv±) (= 5.3 X lO^G cm s~^) with the fiducial temperature. 
(Middle) : The same data are plotted in l//tot — "WiAU plane. Dotted line is the result of Equation J4} adopting the similar conditions to 
those considered in WS91 (see text). (Bottom) : The same data are plotted in -B^,© — fiAU plane. 
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2. FORMULATION 

I derive a simple relation which determines the SW speed near the earth from conditions on the solar surface based 
on a basic energy conservation relation. I consider an open magnetic flux tube measured by heliocentric distance, r, 
which is anchored at the solar surface, r = Rq. Under the steady-state condition, the energy equation becomes 



pv{- 



7 



7-1 



-RT 



GMr. 



(v X B) X B + Fc 



QR = 0, 



(1) 



where p, v, T, B Fc and are density, velocity, temperature, magnetic field strength, conductive flux and radiative 
cooling, respectively. R, 7, G and Mq are respectively gas constant, a ratio of specific heat, the gravitational constant 
and the solar mass. The term involving B denotes Poynting flux under the ideal magnetohydrodynamical approx- 
imation. The cross section of the tube is assumed to expand in pr oportion to r'^f{r), where fir) is a super-radial 
expansion function {/{Rq) — 1 and /(riAu) = ft^t'l ljKopp fc Holzen 11976^1 . Note that divergence of an arbitrary 
vector. A, becomes 'V • A — -^-^(fr^Ar). 

The Poynting flux term can be divided into two parts, —j^{v X B)xB = j^{—BrSB±Sv±+VrSB'j^), where subscript 
r denotes the component along the flux tube and _L indicates the tangential components; SB± and 6v± are amplitudes 
of transverse fluctuations of magnetic field and velocity. The first term indicates shear of magnetic field, corresponding 
to Alfven waves, and the second term denotes advection of magnetic energy. 

Following FSZ, I consider the energy conservation in the flux tube between at the solar surface (r = Rq) and at lAU 
(r = 215i?0). At the surface, besides the gravitational potential, dynamical and magnetical energy associated with the 
surface convection gives an important contribution. Here I rearrange these terms concerning the inputs of energy by 

the convection into two parts, — X B) X B]r + pv{^ + ^^^RT) = — j;^Br5Bi_5v±^ + Fh, namely incompressive 
part (Alfven wave), —j^BrSB±6v±, and compressive part, F^. Note that the magnetic energy term {vrSB'^/An) is 
included in Fh. At lAU the dominant term is the kinetic energy of SW (FSZ). Then, the energy conservation in the 
flux tube gives 



pvr^f tot 



r=lAU 



Br{SBj_Svj_) 

47r 



GMr: 



r=Rn 



lAU 



drr fqji, 



(2) 



where (• • •) denotes time-average. Thermal conduction does not appear explicitly because it only works in redistribution 
of the temperature structure between the two locations. 

The reason of the decomposition of the energy injection terms into the two parts is that their dissipation characters 
are different. Generally, the dissipation of the Alfven wave is slow since it is hardly steepen to shocks. Therefore, 
Alfven waves propagate a long distance to contribute to the heating and acceleration of the SW around ~ a few 
to ^ lOi?0 (SI06). On the other hand, compressive waves an d turbulences denoted b y Fh are more dissip ative so 
that they only contribute to the heating in the chromosphere (|Carlsson fc Steiri1ll992() and the low corona ijSuzuki I 
1200^ . Most of the energy of Fjj which dissipates in the corona is lost by downward thermal c onduction towa rd the 
chromosphere, which finally radiates away in the transition region and the upper chromosphere ijHamme ^ 11981 . The 
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Fig. 2. — Schematic picture of SW in a magnetic flux tube which is super-radially open. Br-,© is proportional to Poynting flux input 
from the surface, /tot determines adiabatic loss in the flux tube. Therefore, the kinetic energy of the SW in the outer region is inferred to 
have positive dependence on Br.Q and negative dependence on /tot. 
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rest (~ 10%) of the energy is transferred to enthalpy flux l|Withbroelll988f) to keep the hot corona with T > lO^K. 
The radiative cooling (last term in Equation I^J is also only efficient in the low corona or below where the density is 
sufficiently high. Therefore, subtraction of the radiative loss from Fa gives 'effective' coronal temperature, Tc, namely 



Finally, we have a conservation equation 



pvr^f tot 



r=lAU 



Br{SB±SvA 

An 



pv 



7- 



-RTc 



(3) 



r=R.r. 



The second term {pvr^ ^^RTq) on the right hand side is evaluated at the base of the corona in the strict sense since 
it implicitly considers the energy balance at the transition region (see §3). However, I use the location, r = i?Q, 
because the distance between the photosphere and the corona is much smaller than Rq. A conceptional novelty of 
the present formulation is that I explicitly include the Alfven wave term which was neglected (FSZ) or parameterized 
in more phenomenological ways (WS91; Schwadron & McComas 2003) in previous works. Note that —5B_\_5va_{> 
for outgoing Alfven waves) is a conserved quantity of Alfven waves which propagate in static media if they do not 
dissipate^. Therefore, —{SB±Sv±) is a measure of dissipation and reflection of Alfven waves in the chromosphere and 
the low corona where the gas is almost static, and the results of numerical simulations (SI05; SI06) can be used for 
the Alfven wave term in a straightforward manner. The physical meaning of Equation Q is clear; the kinetic energy 
of the SW at lAU is determined by positive contributions from the input Alfven wave energy (first term) and the 
thermal pressure of the corona (second term) and a negative contribution due to the gravitational potential well (third 
term). 

Rearranging Equation by using the mass conservation relation, [pvr'^ ftot]r=iA\j = [pvr'^]r=RQ, we can derive a 
more direct form which predicts SW speeds in the vicinity of the earth from the physical conditions on the solar surface 




VlAV = -1/2 X - 



= 300(km/s 



-fXQ 



B 



r,0 



47r(pwr2)iAu /tot 



{SB^SvA 



7- 



1 



RTc~ 



Rr. 



(4) 



-{5B^5vA 



8.3 X 105(cnis-iG) 



/tot 



3.4 



_7_ 
1.1 



0.1 



7-17 V106(K) ) 



4.2, 



where I use the observed mass flux at lAU, (pw)iAU — 5.4 
SWs with different speeds (e.g. Aschwanden, Poland, & Rabin 2001) 
estimated from observed granulation motions at the photosphere as 5vi_ 



10 ^^(g cm ^), which is almost constant even in 



The velo citv amplitude at the surface can be 
~ Ikm/s ijHolweger et al.ll978j) . The magnetic 
amplitude is derived from 5va_ and the photospheric density p ~ lO^^g cm^^ as 5B± = —5v±y/Anp ~ — llOG. Then, 
— {5B^5v\_)q ~ 5.5 X 10®(G cm s~^), where a factor of 1/2 is included due to the time-average. This value is for the 
case in which all the Alfven waves from the surface propagate into the SW region. In the real situation, they suffer 
reflection in the chromosphere. SI05 shows only ~ 15% of the initial energy propagates outwardly to contributes to 
the heating of the coronal and SW plasma. Therefore, we adopt —{5B±5v±)q ~ 8.3 x 10''^(G cm s~^) as a fiducial 
value. As for the thermal pressure, I assume Tc — lO^K as a t ypical coroiial tem perature. 7 should be larger than the 
adiabatic value (7 — 5/3) because of the thermal conduction ijSuess et alJll977j) : I consider 7 = 1.1 in this paper. 

3. RESULTS AND DISCUSSIONS 

I plot relations derived from Equation in the top panel of Figure 1. The fiducial case with {6B±6v±) = 8.3 x 10^(G 
cm s~^) and Tq = lO^(K) (solid line) explains the observed trend quite well. The viau ~ Br^s/ftot relation reflects 
the Alfven waves which accelerate the SW in expanding magnetic flux tubes. S^,© determines energy flux of the 
Alfven waves (oc [BrSB±Sv±\Q). /tot controls 'dilution' of the energy flux; in a flow tube with larger /tot more energy 
is used to expand the tube rather than transfered to the kinetic energy of SW. Thus, the positive dependence on 
Br^Q and the negative dependence on /tot are naturally derived. The result does not depend on different dissipation 
processes of Alfven waves because I only consider the SW speed at the sufhciently distant location (r = lAU) where 
the wave energy are mostly dissipated. 

The top panel of Figure 1 also exhibits that i?r,0//tot is the most important parameter in determining the SW speed 
and that other physical conditions on the solar surface should be similar. Almost all the data are between dot-dashed 
and dashed lines which are the results of cases adopting slightly larger Tc(= 1.5 x lO^K) and smaller {SB±6v±) 
{— 5.3 X lO^G cm s~^) than the fiducial case. Particularly, the difference of the wave amplitudes {5v± cx ^ {SB±6v±)) 
between the solid and dashed lines are only 20%. This indicates that the amplitudes of Alfven waves at the surface 
should be very similar in different coronal holes. The observed data are from not only polar coronal holes but mid- 
latitude and equatorial coronal holes, some of which are located near active regions. Therefore, one can infer that the 
amplitudes could vary a lot since the circumstances are quite different. However, the observation seems to favor the 



Br-^{5B^Svj_), where we have used 



^ This is derived from conservation of the wave energy flux, = -ppj -£:{fr^ BrSB±Sv±) 

conservatio n of magnetic fi eld, Bfr^ =const. Incidently, in moving media wave action should be used as a conserved quantity instead of 
energy flux lJacQues1ll977l) . 
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constancy of the Alfvenic fluctuations at the footpoin ts. Although it is very difficult to observe Alfvenic motions of 
field lines on the solar surface at present ijUlrich Ill99fil ^. this can be observationally studied in the very near future by 
Solar-B satellite which can stably observe fine-scale motions of surface magnetic fields. 

Let me com pare the present analysis with a model calculation for the viau — l//tot correlation by WS91 (see also 
iWang I l|1993')l Although it is not simple to compare both since the assumptions adopted in WS91 are different from 
mine (for example, WS91 fixed the coronal base density, while I adopt the constant mass flux at lAU), I can derive a 
relation of uiau — l//tot from Equation Q by using similar constraints to those considered in WS91. They adopted 
a constant field strength at lAU, Br.iAU = 3 x 10~^G, and assumed a constant energy flux (= 1.5 x lO^erg cm~^s~^) 
at the coronal base. In the middle panel of Figure 1, I present the result with fixed i?r,o//tot — 1-4, corresponding to 
i?r,o = 3 X 10~^G, and energy flux of Alfven waves, Br,Q{SB±Sv±)/4Tr = 4 x lO^erg cm~^s~-^ {i.e. {SB±Sv±) oc f^l) 
in Equation Q by the dotted line. Note that I need the larger energy flux because the inner boundary is not the 
coronal base but the photosphere. The figure shows that the dotted line follows the average trend of the data and 
the result of WS91 is reasonable. However, I would like to address that some data are located away from the main 
viAU — l//tot trend and they can be explained in a unified manner by taking into account i?r,o- 

The relation of Equation se eminglv contradicts to the reported anti-correlation of the coronal temperature and 
the SW speed (|Geiss et all 119951) . This is because my treatment of the thermal processes near the surface is too 
much simplified; the complicated energy balance from the chromosphere to the low coro na is represented only by 
the 'effective' enthalpy, ^i^RTq. The formulation for the temperature- velocity relation bv lSchwadron fc McComasI 

l)2003j) is complementary to the present formulation. In lSchwadron fc McComasI l)2003j) the detailed energy balance 
at the transition region is taken into account, while they assumed a constant input of the Alfven wave energy flux 
which I investigate more in detail. 

In this letter, in order to focus on the SW speed, I simply apply the observed (almost) constant mass flux at lAU 
when deriving the relation for the SW speed. For self-consistent treatments, however, it is important to study how to 
determine the mass flux of the SW not only by numerical simulations (e.g. SI06) but by simple models. 

I think that Equation Q is applicable to SWs during both sunspot minimum and maximum phases because it is 
derived based only on the simple energetics. At present, however, the observed data (Kojima et al.2005) which I use 
for the comparison are only during the sunspot minimum phase (1995-1996). In order to study the generality of the 
derived relation, comparisons with SW data during different phases (Fujiki et al.2006) are important. One should be 
careful that /tot which should be used for the prediction is the actual super-radial expansion factor from Rq to lAU, 
while in most cases, including Kojima et al. (2005), /tot is observationally estimated from the comparison between at 
Rq and at the source surface (2.57?©) based on the potential field-source surface method. Errors due to this method 
could be non- negligible if the potential approximation becomes worse and/or if one separately treats flux tubes in a 
coronal hole (WS90). (In this sense, Kojima et al.2005 as well as I use i?r,o//tot instead of field strength at the source 
surface.) Thus, the precise determination of coronal magnetic fields (e.g. Linker et al.l999) is important for reliable 
forecasts of SW speeds from the relation of Equation Q). 

The author thank the solar wind group in Nagoya-STEL (Kojima, M., Tokumaru, M., Fujiki, K.) for fruitful 
discussions and providing observational data. The author also thank the referee for constructive comments. This work 
is supported by the JSPS Research Fellowship for Young Scientists, grant 4607. 
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